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Abstract. The production of hadrons in relativistic heavy ion collisions is studied 
using a statistical ensemble with thermal and chemical equilibrium. Special attention 
is given to exact conservation laws, i.e. certain charges are treated canonically instead 
of using the usual grand canonical approach. For small systems, the exact conservation 
of baryon number, strangeness and electric charge is to be taken into account. We have 
derived compact, analytical expressions for particle abundances in such ensemble. As 
an application, the change in K/tt ratios in AGS experiments with different interaction 
system sizes is well reproduced. The canonical treatment of three charges becomes 
impractical very quickly with increasing system size. Thus, we draw our attention to 
exact conservation of strangeness, and treat baryon number and electric charge grand 
canonically. We present expressions for particle abundances in such ensemble as well, 
and apply them to reproduce the large variety of particle ratios in GSI SIS 2 A GeV 
Ni-Ni experiments. At the energies considered here, the exact strangeness conservation 
fully accounts for strange particle suppression, and no extra chemical factor is needed. 



1. Introduction 

The statistical models, based on the assumption of the quick natural entropy 
maximization before the freeze-out in the relativistic heavy-ion reaction, have been 
succesful in describing the total hadron abundancies in various experiments, see e.g. 
|U . The usual method of grand canonical ensemble is appropriate in case of very large 
interaction systems. However, the systems studied experimentally are rather small, so 
the exact conservation laws should be taken into account. In such approach, the thermal 
analysis agrees well even with as small systems as induced in CERN e + e~, p — p and 
p — p experiments . 

For larger systems, taking into account the simultaneous conservation of baryon 
number B, strangeness S and electric charge Q is more complicated. We review here the 
method developed, and the application to AGS E802 p— p, p — A and A — A experiments 
in the systematic manner to find if the increase in production of kaons compared to 
pions can be due to finite volume effect incorporated in relativistic canonical statistical 
mechanics @. 

If the net number of any relevant charge becomes large, the task of computing 
abundancies in an ensemble respecting exact conservation of B, S and Q becomes 
difficult. Letting B and Q fluctuate in the grand canonical manner, we may focus 
to strangeness conservation. As an application, we have shown that the particle ratios 
in GSI Ni+Ni experiments at 2 GeV A agree well with the statistical predictions ||. 
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2. The model 



In the usual grand canonical partition function, the states respecting exact quantum 
number conservation are included, but only the conservation of expectation values is 
demanded. By denoting the set of conserved quantum numbers by {C{}, the canonical 
sum of states, Z{ Ci }, can be projected from grand canonical one, Zc |J. In case of N 
U(l) internal symmetries, the projection takes the form 

Z G (T,{A Ci }). (1) 

Here we have assigned a group angle 0j and a Wick -rotated fugacity factor A, = e 1 ^ 
for every exactly conserved charge. 





2-1- Z B s,q 

Let us quote first the expression for the canonical partition function respecting strict 
conservation of B, S and Q. After putting in the set of overall charges, {Ci} = {B, S, Q}, 
and performing some algebraic excercise 0, we find 

Z b ,q,s(T, V) - 

X ^--B+n 2 +n 3 +?i4+?i5+n 6 +n 7 (2A/' n ) 

x ^-Q+ni+n2-n 3 +n5-n 6 +2n 7 (2A/' 7r c ) (2) 
X ^-S , +ni~n4-ri5-n 6 (2A/} S 'o) 

x J m (2^c)/ n2 (27V p )J n3 (2jV A -) 

x J„ 4 (2M A )I n5 (2A/" S+ )/ ne (2A&-) (2A/" A+ + ) • 

Here we have used the modified Bessel functions I n with arguments 2Ai, where Mi is 
the sum of one particle partition functions for particles carrying same set of quantum 
numbers, {Bi, Si, Qi}, as the hadron i. Zq is the partition function for particles with 
vanishing quantum numbers considered. In the expression above, we have omitted the 
contribution from hadrons with \Si\ > 1, but the generalization is straightforward 0. 
The mean abundance of hadron i in the system is 

^jy.^ B—Bi,Q—Qi,S—Si ^7! 

Zb,q,s 

The evaluation of the canonical partition function with three simultaneously 
conserved quantum numbers becomes numerically very time consuming for large values 
of B. So far, for systems with B > 20 we have been forced to resort to the grand 
canonical treatment. 
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2.2. Z s 

In large systems, such as Ni+Ni, the grand canonical treatment of baryon number 
and electric charge is justified. However, the exact strangeness conservation is still 
mandatory in theoretical considerations. If we only include the particles with \Si\ < 1, 
and require vanishing net strangeness, we end up with ||, [TJ 



Z S=Q (T, V, X B , A Q ) = ZM2yJNiM-i), (4) 

where the iV^ are the sums of grand canonical one-species partition functions of particles 
carrying strangeness Si defined by 

Z l = X b X Q^ jf dpp^V^. (5) 

The equation (f|) can be generalized to include the particles with higher strangeness 
content as well, see || . Now the mean number of hadrons i is 

The nonlinear volume dependence of strange particle production rates is visible 
in the coefficient in equation (|5]). Whereas in grand canonical formalism the 

abundance of kaons is linear in volume (Zl), the Ii/Iq term gives an additional coefficient 
of first order in volume in small system limit. This nonlinearity decreases smoothly with 
increasing volume, until in thermodynamic limit it vanishes. 

3. Applications and discussion 

3.1. K/tt ratios in AGS experiment E802 

We have applied the Zb,s,q in comparison with experimental K/tt production ratios 
reported by E802 collaboration, see table |I|. In figure [l] we find the theoretical curves 



Table 1. Experimental results reported by the E802 collaboration. B and Q are 
calculated using geometrical considerations, see || 



Collision K+/tt+ 



P - 
P - 
P - 
P - 
Si 



Be 

Al 

Cu 

Au 

Au 



Ref. 



Ref. 



B 



7.8±0.4% 
9.9±0.5% 
10.8±0.6% 
12.5±0.6% 
18.2±0.9% 
19.2±3% 



2.0±0.2% 



2.8±0.3% 
3.2±0.3% 
3.6±0.8% 



3.9 
5.4 
6.9 
9.7 
102.7 



Q 



2.3 
3.1 
3.7 
4.5 
44.0 
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for K + /tt + and K~/7r~ ratios in constant baryon density as functions of baryon number 
to be consistent with experimental results. This shows that the increase in K/n ratios 
along with increasing system size can also be explained without any phase transition or 
the in-medium effect in kaon masses [111], O, [T3] . 



3.2. Production ratios in GSI Ni+Ni experiments at 2 GeV A 

As an application of Z$ the thermal model analysis is compared to experimental hadronic 
ratios in GSI Ni+Ni 2 GeV A experiments. In this case, the widths of resonances affect 
substantially the results, so we have applied the relativistic Breit-Wigner resonance 
shape in computing the phase space integrals (|5|). 

In Ni+Ni system, the isospin asymmetry = 1.04) has to be taken into account 
by introducing chemical potential for electric charge. This parameter, however, is 
eliminated by the simple binding condition between baryon- and charge densities, ub 
and riQ\ 

n B = 2(^)rkj. (7) 

The comparison between model predictions and experimental data is collected in 
table ^. Only the ratio r]/iT does not fit reasonably in thermal model scheme, but 



Table 2. Particle 
experimental results, 
potential is used. 



ratios resulting from Zg 
The best fit value, [is - 



compared to GSI Ni+Ni 2 GeV A 
= 0.72 GeV, for the baryon chemical 



Ratio 


Model 






Data 




R [fin] 


4.2 


3 








T [MeV] 


65 75 


65 


75 


ratio 


ref. 



K+/K" 

K+/7T+ 

<t>/K- 

77/7T 
7T+/p 
7T°/B 
d/ P 



25.7 

0.0071 

0.103 

0.893 

0.008 

0.225 

0.104 

0.129 



22.4 

0.0339 

0.082 

0.895 

0.015 

0.247 

0.108 

0.188 



23.9 

0.0027 

0.276 

0.894 

0.008 

0.224 

0.104 

0.129 



21.1 

0.0132 

0.212 

0.898 

0.015 

0.246 

0.107 

0.188 



21 ±9 

0.0074 ±0.0021 

0.1 

0.89 

0.037 ±0.002 
0.195 ±0.020 
0.125 ±0.007 
0.26 
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it serves as a clear indication of incompleteness of the thermalized ideal gas model 
considered. Apart from that, the other measured hadronic ratios are well described 
by the model. As a clarification, the equal value curves for ratios in (T, plane in 
constant volume are shown in figure [2[ 
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Figure captions 



Figure 1. Thermal model expectations for the production ratios K + /tt + and 
K~ /ir~ at a temperature of 100 MeV and a baryon density of 0.04 fm~ 3 compared to 
experimental results from the Brookhaven AGS. The experimental ratios from Si — Au 
collisions (B ~ 103) is moved to B = 21 for the sake of convenience. 



Figure 2. Curves in the ([ib,T) plane corresponding to the GSI Ni+Ni 2 A GeV 
hadronic ratios indicated. The interaction volume corresponds to the radius of 4.2 fm, 
and the isospin asymmetry is B/2Q = 1.04. 
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